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CN ■ Abstract 



O ! In this article, we are interested in the behaviour of a single ferromagnetic 

^^ ' mono-domain particle submitted to an external field with a stochastic perturbation. 

This model is the first step toward the mathematical understanding of thermal effects on 

a fcrromagnct. In a first part, wc present the stochastic model and prove that the 

associated stochastic differential equation is well defined. The second part is dedicated 

to the study of the long time behaviour of the magnetic moment and in the third part 

we prove that the stochastic perturbation induces a non reversibility phenomenon. Last, 

we illustrate these results through numerical simulations of our stochastic model. 

,S^ I The main results presented in this article are on the one hand the rate of convergence of 

j^ ' the magnetization toward the unique stable equilibrium of the deterministic model and 

on the other hand a sharp estimate of the hysteresis phenomenon induced by the 

stochastic perturbation (remember that with no perturbation, the magnetic moment 

remains constant). 

|J^ I Keywords: convergence rate, stochastic dynamical systems, long time study, magnetism, 

r — ■ hysteresis. 
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O ■ 1 Introduction 



Thermal effects in ferromagnetic materials are essential in order to understand their 
behaviour at ambient temperature or, more critically, in electronic devices where the joule 

k> ' effect induces high heat fluxes. This effect is commonly modeled by the introduction of a 

5h ■ noise at micro-scopic scale on the magnetic moment direction and at the meso-scopic scale 

- - -' by a transition of behaviour. In ferromagnetic materials, the transition between the 

non-linear behaviour and the linear behaviour is managed by the struggle between the 
Heisenberg interaction and this disorder induced by the heating. This model explains the 
critical temperatures such as the Curie temperature for ferromagnetic materials. In this 
context, it is essential to understand the impact of introducing stochastic perturbat ions in 



deterministic models of ferromagnetic materials such as the micromagnetism (see iBrown 
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(|l962l . ll963l ^). 

The understanding of this phenomena is a key point in order to simulate realistic 
ferromagnetic devices such as micro electronic circuits. Furthermore, heating has a real 
effect on the microstructure dynamics in magnets; then, efficiently controlled, the dynamic 
of microstructures could accelerate processes such as the magnetization switching, which is 
the basics of magnetic recording techniques. 



Duri ng the last de c ade, s e veral studies have b een initialized i n seve r al articles by physicist s 
fe.g. iMercer et all (I2OIIII: Izherig et al.l (I2OO3I]: iRaikher et all toOJ) : lAtkinson et all (I2OO3I ): 
Raikher and Steoanovl \2mit ]: IScholz et al.l (|2nnih : iMartinez et al.l (|2nn7l ^: ISmit'hr(|2nnih ). 
but no mathematical models justifying this kind of effects have been developed yet. In this 
article, our goal is to improve the understanding of thermal effects in ferromagnets. To 
achieve this goal, we focus this first study on the dynamic of a single magnetic moment 
submitted to a stochastic perturbation. Our main aim is to characterize precisely the 
dynamic of the moment, giving estimations and general behaviour in long time. In 
particular, we will exhibit an hysteresis behaviour of the magnetization in our model. 



The model we are studying mimics the behaviour of a single magnetic moment /i(t) (function 
from M into 5(]R'^) = {u £ M'^; |/i| = 1}) submitted to an external field b. The dynamic of 
such a system is, at the micro-scale, described by the Larmor precession equation 

m 

Nevertheless, this equation is non dissipative and, in order to make the theoretical study 
easier, we introduce a dissipative part using the Landau-Lifchitz equation 

dfj, 
Ht 

where a is a positive real constant and we set the initial condition /i(0) = fiQ £ 5(R^). We 
point out two major properties of this system 

i. ViGM, |^(t)| = 1, 



-/.f A b. 



-/.f A 6 — af-i A {fi Ab), 



a. \/t G 



|(Mt)-&)>o. 



The first property, which is definitely essential, will have to be preserved by the stochastic 
system and the second property is the energy decreasing induced by the introduction of the 
dissipation term. The dynamic of this deterministic system is classical; in fact, one knows 

that lim fxit) = —rb provided that ^(0) 7^ —b. In this work, we will develop such a result 

t-S'Oo \b\ 

for the stochastic system. In order to build this stochastic system, the first question is how 
to introduce the stochastic perturbation in the deterministic system. We want to model 
the thermal effects which are external perturbations of the magnetic moment. In fact, this 
perturbation could be modeled has an external perturbation field. In the sequel, we will 
choose to build a stochastic system by perturbing the external field with a Brownian motion. 
We will write down 



-fit A{bdt + £ dWt) 



a^H A utAib dt + e dWt) 



where e is a strictly positive real number and W a standard Brownian motion with values in 
M^. But, an easy computation of d(|ltp) using Ito's formula shows that the process Y will 
not stay in S'(IR^), then, in order to preserve this essential behaviour, we have to renormalise 
the previous equation and set 

Yt 

Given this system, we prove the following results 
i. fit • b > \b\, a.s., 

ii. limt s>oo Vt IE[||^| — fJ-t • b\] exists and we can compute its value. 

Note that these results are only valid for a > 0, as for a = 0,it is easy to show that the function 
e(t) = E(//t • b) satisfies the following ordinary differential equation e'{t) = — ^(t)-^^. Hence, 

e{t) = ft4 > 0. This contradicts the a.s. convergence of /Uj to t^. 

in. When n is submitted to a time varying external field, an hysteresis phenomenon appears. 
If we consider b G 5(M^) and let b linearly vary between +6 and —b over the time interval 
[0, r], then E(;U( • b) is bounded from below by ,}_ for t < T/2 where c is a constant 
depending only on e and a. 

First, we make precise the derivation of the stochastic model. Then, we lead a detailed study 
of its asymptotic behaviour and in particular we point out an hysteresis phenomenon. This 
phenomenon is obtained by slow variations of the external field such that the dynamic of 
relaxation of the magnetization toward this field becomes instantaneous when the speed ratio 
of the external excitation goes to zero. The results shown in this article are finally illustrated 
by numerical simulations. 

2 Model and notations 

Let {Q, A, P) be a probability space. We consider a standard Brownian motion W defined on 
this space with values in M^ and denote by F = {J-t)t>o its natural filtration augmented with 
Pnull sets. 

Let 5 S M^ be the magnetic field. We model the 5(M'^)-valued magnetic moment process 
fi = (//t)t>o by the following coupled stochastic differential equation (SDE in short) 

-/if A{b dt + e dWt) - afit A fit A {b dt + e dWt) 

^ (2.1) 



Wt. 
yes 

where a > is the magnitude of the damping term and e > is the magnitude of the noise 

term. 

The term fit A dWt in (|2.1|) is naturally defined by introducing the antisymmetric operator 

L : W^ I — > M'^^^ associated to the vector product in R^ 

—x^ x^ 
L(x) = \ x^ —x^ 
-x^ x^ 



Hence, for a 3-diinensfonal F— adapted process H satisfying /q \Hs\ ds < oo a.s. for all t > 0, 
the process /q Hg A dWg is defined by /q L{Hs)dWs, which is a standard multi-dimensional 
Ito stochastic integral. When dealing with the differential expression, we will either write 
HtAdWt or L{Ht)dWt. 

Notations: 

• For a and b in R^ we denote by a • 6 their scalar product, a ■ b = J2i=i ^*^*- 

• For a in M'^, we denote by \a\ = ^/a ■ a the Euclidean norm of a. 

• We like to encode elements of M'^ as column vectors. For x G M.^, x* is a row vector. 
Similarly, we use the star notation "*" to denote the transpose of matrices. 

• li H = {Ht)t>o is a 3-dimensional F— adapted process satisfying /q \Hu\'^du < oo a.s. 
for all t, we may write /q Hu ■ dWu for X]i=i /o ^u dW^ and use the differential form 
Ht • dWt for J2f=i Hi dWl 

3 Main results: long time behaviour 

3.1 First properties of the magnetic moment /i 
Proposition 1. Let iY^^) be a pair of processes satisfying (|2.1|) . then 

d\Yt\^ = 2e^{a^ + l)dt 



and therefore \Yt\ = \j2e^{pP' + l)i + 1 is non random. 

Remark 2. The fact that \Yt\ is non random, is definitely essential in all the following com- 
putations. In particular, we deduce from this result that \Yt\ has finite variation. 

Proof. Using Ito's lemma we have 

3 3 






d\Yt\^ = 2Yt -dYt + Y, d{Y\ Y')t = ^ d{Y\ Y')t, 



where we have used the fact that Yt and dYt are orthogonal. But, using the identity aA{bAc) 
(a • c)b — (a • b)c, we have 

dYt = e[- (nt A dWt) - aiifit ■ dWt)i^t - (f^t ■ f^t)dWt)] +...dt 
= eA{fit) dWt + ... dt, 

where we have set A{n) = al — a{fifj,*) — L{n) and used \iit\ = 1- Thus, 

d{Y,Y)t = e^AA*{fit)dt 

a^I - a'^ifitfJ't.) + aL{fit) - a^{^tlh) + "^(/^tA*t )(A*t/^t ) - »ilJ'tlJ't)L{tJ-t) 



£2 



- aL{fit) + L{nt)if^tlJ't) - Hl^t)Lifit) 



dt 



where we have used L{n)fi = 0, L*{fi) = —L{fi) and again /-^/Ui = 1. Thus, for each I < i < 3 
we have 

d{Y\Y% = e\a^ - a\,I,f + ^(L,fc(/ii))2] . 

fc=i 
Then, summing over i, we get 

3 

d\Yt\^ = e2[3a2 - a^ \^it\^ + ^ (L,,(//i))2] dt 

= e^[Za^ - a^ |/it|^ + 2 |^t|^]di, 
The result ensues by remembering that |/ii| = 1. ■ 

With the help of Proposition [H we can establish the SDE satisfied by the one dimensional 
process (/Uj • b)t- We introduce the function 



/i(t) = |y(t)| = ^2e2(a2 + l)t + l. (3.1) 

Since |^| is non random, we deduce from Equation ()2.ip that 

= - ^w \ dt - —— {fit /\{bdt + e dWt) + a/i* A ^it ^ {b dt + e dWt)) 

, Mt/t'(t) + fit/\b + ajutif^t ■b)-b) ^ (Ti \ , I * t\\m;(^ ^q o^ 

dm = — dt - — {L{nt) + a{^It^it - I)) dWt (3.2) 

By taking the scalar product with 6, we get 

d{^,t . b) = -(/., • b)^dt - ^ ((^, • bf - 16|2) dt 

/i(t) 

d{^,t . b) = -(/., • b)^dt - ^ ((^, • bf - 16|2) dt 

We may call this equation the SDE satisfied hy ^ ■ b whereas it is not an SDE properly 
speaking since the r.h.s member actually depends on all the components of /i^ and not only 
on /if • b. Nonetheless, we may use this abuse of terminology throughout the paper. 

Remark 3 (Remark on the existence and uniqueness of solutions to Equation (j2.1|) ). Let us 
consider the following coupled SDE 

dYt = -fit Aibdt + e dWt) - a^t A fit /\ {b dt + e dWt) (3.4a) 

dut = ^^T^T dt - — {L{fit) + aif^tfJ-t - I)) dWt (3.4b) 

Yo = fio G SiW") 



This system is actually decoupled as the SDE on ^ is autonomous (this has only been possible 
because \Yt\ is non random). The existence and uniqueness of a solution to Equation (|3.4p 
boil down to the ones of Equation ()3.4bp . By computing fi(|/xjp), we deduce that if there exists 
a solution (i to Equation ()3.4bp . |/it| = 1 a.s. for all t. Hence, it is sufficient to check the 
standard global Lipschitz behaviour of the coefficients on M_(_ x 5(M^) to prove existence and 
uniqueness of a strong solution to Equation ()3.4bp . 

We have already seen above that if a pair {Y, fi) is solution of Equation i\2.1\] . it also solves 
Equation (|3.4p . 

Conversely, if {Y, fx) is the unique strong solution of Equation ()3.4p . it is clear that |lj| = h{t) 
by following the proof of Proposition [J\ and moreover the computation of d{Yt/\Yt\) shows 
that the process {Yt/\Yt\)t solves the same SDE as n, hence for all t fit = Ty] ^■^- This 
last argument proves that Equations (|2.ip and (|3.4p have the same solutions. Therefore, we 
deduce that Equation ()2.ip admits a unique strong solution denoted by (1", /i) in the sequel. 

3.2 Almost sure convergence 

In this part, we prove the almost sure convergence of /ij to b/\h\ when t goes to infinity. This 
is achieved by studying the pathwise behaviour of the process (/i^ ■b)t. 

Theorem 4. lim ^uj • 6 = |6| a.s. 

t — >oo 

To prove this result, we need a preliminary result stating that the stochastic integral in 
Equation ()3.4bp actually vanishes at infinity. 

Lemma 5. 

sup/ -TT^{~ l^u^b + a{{^u-h)iiu-b)) ■ dWu<oo a.s. 
t Jo h(u) 

Proof of Theorem \^ From Lemma El we know that 

sup / TT^{ - fj,u/\b + a{{fj,u ■ b)fj,u - b)) ■ dWu < oo a.s. 
t Jo h(u) 

Hence, we can define for all t > 

/•oo s 

Xt = Hfb- / — - ( - /i„ A 6 + «((//„ • b)iJ,u - b)) ■ dWu 

Let |6| > 5 > be chosen close to 0, there exists T such that for all t > T, \Xt — fit ' b\ < S. 
Moreover from Equation (|3.3p . we can deduce that for all t > s > T 

Xt-Xs= f -flu ■ b^ + T^iifiu ■ bf - \b\^)du. (3.5) 

Js h(u) n(u) 

Let r/ < |6| be chosen close to |6|. On the set {0 < /x^ • & < r/} we have \fiu A 6| > \b\ — rj^. 
Thus, 

,,/i'(n) a , ,,2 ,,,/i'(w) \b\^-ri^ 

-{^^u■b)-^ + ---\^iu^b\^>-\b\-^ + a^-^ 



h{u) h{u) h{u) h(u) 



We can always choose T such that for all u > T, 



^nh'iu) \b\ 

— \b\ , , , +Q 



h{u) 



h{u) 



V\ \b\ 
> a 



2h{u) 



On the set {/i„ • 6 < 0} we have 
, ^,/i'(u) 



a 



h{u) h{u) 



W^b\'^>{-{^Mu■h) + a\^lu^b\'^)- ^'^' 



h{u) 



> min(|5|,a|6n-- ^ ^ 



2 h{u) 



The last inequality conies from the fact that if 7r/2 < x < 37r/2, we have either — cos(rE) > 

^/2 or |sin(x)| > i/2/2. 

Therefore, there exists T > T, such that for all u> U, on the event fi^'b < rj, —{fiu ■ 6) r/ ■: + 



>'(«) 



h(u) 



■j^ |^„ A 6| is bounded from below by c-^jw^ where c is a positive real constant depending 
on T]. Therefore, we deduce from Equation 



Xt - X, > 



> 



1^1 
1^1 



*/i'l 



u 



h{u) 
* h'{u) 



'^{x^>r,-S}du + c 



h'{u) 
h{u) 
* h'{u) 



^{i^u-b<v}du 



{Xn<ri^5} 



du 



ft h'(u) 



h{u) -y^^>^-°i"'- -j^ h{u) 

While choosing 5 and r/, we can always ensure that Q < r] — 25 < \b\. 

Assume that Xg < tj — 26. Then, there exists t for which the therm J^ "uul ^{Xu<ri-5}d'u has 
driven X back into the cap such that Xt > rj — 26. Hence, we can assume that Xs > rj — 26. 
Either, for all t, Xt > rj — 26, or thanks to the continuity of X there exists to for which 
XtQ = rj — 26 and within a finite time the second integral drives X back into the polar cap 
such that at a some time Xt > rj — 6 and for all u & [to, t] rj — 6 < X^ > r] — 26. This reasoning 
enables to prove to Xt > r] — 26, ie. Ht ■ b > r] — 3(5. As r/ can chosen arbitrarily close to |6| 
and 6 arbitrarily small, this proves the almost convergence of /Lij • 6 to |6|. ■ 

Proof of Lemma From Doob's inequality we have 

r-t 1 



E 



sup 

t 



K'^) 



{- IJ,uAb + a{{fiu ■ b)nu - b)) ■ dWu 



< 



1 



h{uy 



:E 



fj.u-b\ ) (1 + a^)du. 



(3.6) 



Now, we will prove that the r.h.s is almost surely finite. 

From Equation ()3.3p . we get after integrating and taking the expectation for all f > 



E[/ii • b] - E[/io • b] 
¥.[jifb]'h{t)+¥.[iifb]h'{t) 

E[//*-6]-^E[/.o-6] 



a 



n^^u-b]-^-—- 
h(u) h(u) 



E 



{^Jiu■bf-\bY 



du 



E 



a 

W) 

a 



\b\' - if^t ■ bf 



E 



\b\'-{f,u-bf 



du 



Hence, we deduce that 



t h{t) Jo 



E 



\b\'-{f,u-b)' 



du < —^ = K. 



a 



Let us consider the upper-bound in Equation (|3.6|) truncated to t and perform an integration 
by parts to obtain 



1 



/i(n)2 
1 



E 



{\bf-W-bf) 



du 



h{uf Jo 



E 



i\b\'-\f^.-b\') 



dv 



1 /•* h'(u) , 

< K—- + 2k I -^--TTrdu 



1 



h{t) Jo h{u) 

/I /I 

< K [ ^T-r + 2 ' 



< 



+ 



2k 



t Oh' 



2h'{u) r 



/i(n)3 Jo 



E 



m'-\t^vb\') 



dvdu 



\h{t) \h{0) h{t)JJ - h{0) 
This proves that the r.h.s of Equation (|3.6|) is finite and ends the proof of Lemma O ■ 

3.3 Convergence rate 

In this section, we are in the behaviour of h{t){\b\ — fifb). In particular, we estabhsh the rate 
of decrease of the L^ norm of \b\ — fj,t • b to zero. 



Theorem 6. hm E[h{t) \\b\ - m ■ b\] 

t — >oo 



8^(1 + a^) 
2a 



Proof. As \b\ — fit ■ b > 0, the L^ norm boils down to a basic expectation. 
Let us define ^t = \b\ — fJ-t • b for t > 0. From Equation ()3.3|) . we get 



. h'it) , ,,, h'(t) , a 



h{t) 
h'it) 



h{t) h{t) 



{{f^t ■ bf - \b\' 



dt 



hit) 



fit/\b + a{{nt ■ b)fj,t - b)) ■ dWt 



dt + \h\^dt 



-Um-it)dt 



■{-fit/\b + a{{ixt ■ b)nt - b)) ■ dWt 



^'h{t) ' ' h{t) h{ty ' ' ^'' h{t) 

If we introduce Zt = h{t){\b\ — fit ■ b), we can write 

dZt = {h'{t) \b\ - q6(2 |6| - 6)) dt + e{-fitAb + aiifit ■ b)nt - b)) ■ dWt 
From the dynamics of Y, we deduce that E[Zt] solves the following differential equation 

E[Zt]' = \b\h'it)-a{2\b\E[^t]-m!]) 
E[Zt]' = \b\h'{t)-a2\b\^ + aE[C^] 



E[Zt]' = \b\ h'{t) - a 



2\b\ 



e2(a2 + 1) 



h'{t)E[Zt]+aE[il 



ft 2a\b\ ,11 ^ , \ ' 



(|6|/i'(t) + aE[e?])e-^o;^f^'^'W'^- 



Now, we can integrate the previous equation to obtain 

Jo 

+ a 






+ ae ^^('^'+1) / lE[C«]e" (" +1' du 

Jo 

From TheoremlU we know that ^t tends to a.s., therefore the bounded convergence theorem 
yields that hm„ ^^o ^[Cu] — 0- Hence, as h{t) tends to infinity with t, it is easy to show that 

t S-OO In 



Then, we can deduce that 



t — ^oo la 



As a corollary of Theorem El we can prove the following results using Markov's inequality. 
Corollary 7. For all < P < 1/2 and ?? > 0, P(t^(|6| - Ht ■ b) > r]) — ^ 0. 

4 Hysteresis phenomena 

In this section, we want to study the impact of the stochastic perturbation on the reversibil- 
ity of the system; we are wondering whether the stochastic part may induce an hysteresis 
phenomenon. In order to observe this, the particle is submitted to an external field linearly 
varying from +6 to —b where b G 5(]R^) and with constant direction and bounded modulus. 
We have seen in Section [3] that when the external field is fixed, the magnetic moment /U 
asymptotically stabilizes along this field. If the external field varies sufficiently slowly com- 
pared to the stabilization rate of /z, we expect that // will take different back an forth paths 
when the external field switches from +b to —b and then from —b to -|-6: this is the non 
reversibility property of the system. 

In order to highlight this property, we will study the evolution of a suitably rescaled system 
on the time interval [0, 1] and show that the average back and forth paths of ^t • 6 can not 
cross at the point t = 1/2. 

We consider a two time scale model: a slower scale for the variations of the external field and 
a faster scale for the Landau Lifshitz evolution of the magnetic moment. 
Let ry > be a fixed time scale and b G 5(]R'^) the direction of the external field. We define 
the external filed W linearly varying between +b and —b on the interval [0, l/ij] by 

5'?(t) = {l-2tr])b fort G [0, 1/ry] 



We assume that the magnetic moment /x'' is affected by b^{t) according to the following 
equation for t E [0, 1/77] 

' dYt'' = -lil A {U^ (t) dt + edWt)- a/x? ^ /^t ^ {^"^ (*) dt + e dWt) 

JS = b 

In order to work on the interval [0, 1], we introduce rescaled versions of both the external field 
and the magnetic moment defined for t S [0, 1]. 

Using the time scale property of the stochastic integral, we can write 

dZt = -At A (b{t) -dt + e dWt/rA - aXt A At A ib{t) -dt + e dWt/r, 

From the scaling property of the Brownian motion, we know that (y^VFt /^) is still a Brownian 
motion. So we get 

dZt = -At A I b(t) -dt + e dWt | - aAt A At A | bit) -dt + e dWt | 

\ It] y/n J \ V y/n J 

It is important to notice that the factor r] acts as a time scale parameter for the deterministic 
part, but that the corresponding scaling parameter for the stochastic part is ^/rj. 
Following the proof of Proposition [U it is obvious to show that d\Zt\ = 2(1 + a'^)e'^/r] dt. 
Then, we introduce for t £ [0, 1] 



The Ito formula applied to (At • b)t yields 

h^'(t) 1 — 2t 

diXt ■ b) = -(At • b)-—j-^dt + a—-—- lAt A 51^ dt 



^ ((At A dWt) ■ b + a{\t ■ b){Xt ■ dWt) - a(b ■ dWt)) 



Vvh"{t) 



The stochastic part vanishes through when taking expectation as in the previous section to 
find 

E(At ■b)-Xo-b= -E{Xu ■ b)—^ + a— -— E |A„ A b\^ du (4.1) 

Jo M[u) ■qh^iu) 

Proposition 8. For all t £ [0, 1/2], 

E(At • b) > ——- > ., =. 



10 



Proof. Since the process A is pathwise continuous and bounded, it is easy to show that 
the deterministic function e{t) : t i — > IE(At • b) is of class C^. Hence, we can differentiate 
Equation (fO]) 

/. ^ ^.h'Tit) l-2t ,, ,,2 

e'(t) = -e t — 4t- + Q , . ^ lE Af A far 

As i < 1/2, the second term on the r.h.s is non-negative, hence 



e'(t) > -e(i) 



/i'?(t) 



From this inequahty, we deduce that {e{t)h^{t))' > 0, which leads to the following lower 
bound 

e{t) > -r-TT for t < 1/2 



5 Numerical experiments 

In this section, we want to illustrate the theoretical results obtained in Sections [3] and H] using 
some numerical simulations. 

Long time behaviour. We consider the stochastic system (j2.1|) and discretize it with the 

help of an Euler scheme (y, /i), on a time grid with step size 5t > 0. 

Figure [T] shows the long time behaviour of {fj,t ■ b)t>o, for one path of the scheme {Y, fl), with 

time step size 5t = 0.01 and for different values of the damping parameter a. The parameter 

£ is fixed to 0.1, we have taken |6| = 1, and set /xq = —b. The almost sure convergence of fit -b 

to \b\, as stated by Theorem HI is well illustrated by Figure [T] and one can also see how the 

parameter a impacts the characteristic time of the system, ie. the time needed to stabilize 

around the limit. 

Now, we wish to compare the rates of convergence studied in Subsection 13.31 to numerical 

observations. From Theorem El — ,"^ \/t E(|6| — fifb) converges to 1 when t goes to 

infinity. This is illustrated by Figure [2] for different values of a. This figure confirms that 
decreasing the parameter a leads to a decrease of the convergence rate of ]E(|6| — fifb)- 

Hysteresis phenomena. On Figure El we can observe a typical pathwise hysteresis phe- 
nomenon, which not only illustrates Proposition [5] but also suggests that the result of this 
Proposition could well improved by proving a almost sure lower bound (probably for suffi- 
ciently small values rj) . On Figure El the forward path (red curve) is almost stuck to the 
value 1 on the interval [0,1], we could then be tempted to think that the lower bound of 
Proposition [5| lacks some accuracy. 

On the contrary, when r] becomes small, which corresponds to a slower scale for the variations 
of the external field. Figures Sj and \5\ show that the lower bound l/h^{t) becomes nearly 
optimal for t lower but close to 1/2. 
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Figure 1: Almost sure convergence of /xj • 5 with fiQ = —b, |6| = 1, e = 0.1. 

6 Conclusion 

In this article, we analyzed the long time behaviour of a non linear SDE modeling the evolution 
of a magnet submitted to a perturbated external field. The rate of convergence of the mag- 
netic moment is particularly interesting and holds for any dissipation coefficient a > 0. This 
result has been obtained by combining the ODE technique with Ito's formula. The second 
result concerns the hysteresis behaviour of the system induced by the stochastic perturba- 
tion. The combination of these two results illustrated the dissipative effects of the stochastic 
perturbation on the system by giving a first glimpse on thermal effects on a ferromagnet can 
be modeled in the framework of micromagnetism. 
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Figure 2: Convergence of ^"^, Vt E(|6| - fifb) with /xq = -b, 161 = 1 and e = 0.1. The 

horizontal dashed hne is at level one. The expectation is computed using a Monte-Carlo 
method with 100 samples. 
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Figure 3: Pathwise hysteresis phenomena with a = 1, e = 0.005 and r] = 0.01. The red curve 
is the forward path whereas the blue curve is the backward path. 
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Figure 4: Pathwise hysteresis phenomena with a = 1, e = 0.01 and rj = 3.1E — 5. The blue 
curve is the evolution of fit ■ b and the green curve is l/h^{t). 
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Figure 5: Zoom of Figure U] around t = 1/2. The blue curve is the evolution of fit ■ b and the 
green curve is l/h^{t). 
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